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Abstract 

The scaling hypothesis for the coupled chaotic map lattices (CML) is formulated. Scaling properties of the 
CML in the regime of extensive chaos observed numerically before is justified analytically. The asymptotic 
Liapunov exponents spectrum for coupled piece-wise linear chaotic maps is found. 



1 Introduction 

1.1 The Liapunov exponents spectrum 

Partial differential equations describing continuous models and real systems can be discretized 
into a system of coupled maps called as coupled map lattices (CML). In the last fifteen years, a 
close attention has also been drawn to these models in virtue of studies of the generic properties 
of spatiotemporal chaos, Nevertheless, by this time, just a few rigorous analytical results 

on the chaotic dynamics observed in CML have been achieved, 0-0. 



Chaotic behavior in extended systems is typically incoherent, [|12|]. The problem of stability 
in lattice dynamical systems has been extensively discussed in J7|. The long-time limit along 
orbit of the rate of separation of points in the tangent motion is given by the Liapunov exponents 
spectrum ||, 

1 n— 1 

X k = lim -V log | (1) 

t=0 

where (DF f ) is the Jacobian matrix evaluated on the orbit F l . The notation | . . . \ k denotes 
the k-th eigenvalue of the matrix {DF), ordered in such a way that \(DF)\ k _ 1 > \(DF)\ k > 
\(DF)\ k+1 . The spectrum (|I|) characterizes the stability of motion in the phase space. If all 
exponents are negative the trajectory has an attractive fixed point. If the maximal value of 
the spectrum Ao = 0, the attractor is a limit cycle, whereas positive exponents correspond to 
diverging orbits. 

A simple example of CML is given by the 1-dimensional chain 

F(c t x ) = c t x +1 = (l-2s)f(c t x )+e[f(c t x _ 1 ) + f(c t x+1 )], x E Z, t E Z + , (2) 
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of loosely coupled (e < 1/4) piecewise linear maps 

f(c t x +l )=rc t x mod 1, (3) 

which is known to induce the chaotic behavior as the value of the control parameter becomes 
r > 1. Space-time chaos properties in chains of weakly interacting hyperbolic maps has been 
considered in ||. 

The Liapunov exponent of the solitary uncoupled map ([3D is known exactly as a function of 
the control parameter r, Ao = logr. Furthermore, since the correspondent Jacobian in (0) is a 
N x N tridiagonal matrix, its eigenvalues can be computed explicitly for any finite number of 
sites iV provided the boundary conditions (c_i and c/v+i) are known, Hi~0| . For example, chosen 
the zero valued boundary conditions, c_i = c^+\ = 0, one obtains the following Liapunov 
exponents spectrum, 

/ irk \ 

X k = A + log( l-2e + 2ecos^^J , k = l...N. (4) 

For any finite N, the spectrum (Q) is bounded from above by the maximal value coinciding with 
the Liapunov exponent of uncoupled maps Ao- Either of the spectral values Xk < Ao depend 
both on the boundary conditions chosen and the actual value of the coupling strength parameter 
e. In the extended limit, iV — > oo, the boundary conditions can be assumed as an asymptote 
limjv^oo c±n = 0, and the spectrum (|j) turns to a band. In actual numerical simulations of the 
extended limit, the boundary conditions are usually taken as periodic. 

More realistic models for partial differential equations refer to more complicated coupled 
maps lattices for which the Liapunov exponents spectra have not been found. Therefore, one 
is interested in methods which can yield at least qualitative information on them. 

The Liapunov exponents spectra of CML in the extended limit, N — > oo, is of special interest, 



I2|j . Intuitively, if the system is large enough, subsystems that are sufficiently separated in space 



become decorrelated and so contribute to the fractal dimension in proportion to their number. 



This phenomenon usually called as an extensive chaos, [|TTJ, is numerically observed for lattice 
sizes larger than some threshold value depending on the model, ||. 

For large lattices, the number of positive Liapunov exponents Xk > scales like the lattice 
size 

N Xk>0 ~ N/U, (5) 

where N is the volume of a large albeit finite piece L of the lattice, and U is said to be the 
volume of a typical statistically independent block. In accordance with numerical simulations 
data, U depends substantially on the particular model while the relation fl5|) is believed to hold 
generally. The result fl5|) has been observed numerically for some simple coupled map lattices 
(see Jl3| and references therein for a review). Nevertheless, up to our knowledge, the scaling @ 
has not been justified analytically, in particular, the meaning of the statistically independent 
block U is still vague. 



Another fascinating phenomenon observed in the chaotic coupled map lattices ||14j|-|12l is a 



non-trivial collective behavior ("hidden coherence", |JL5||), where there is apparently a unique 
nice invariant measure, PBj, but nevertheless non-uniqueness of nice measures on the space- 



time configurations, because a two-cycle or more complicated dynamic behavior occurs in time. 
The time periodic behavior of measures has been discussed in Jl6| ("asymptotic periodicity"). 



The chaotic dynamics in CML lacks of mathematically rigorous results. Nevertheless, some 
"physical" approaches borrowed from statistical physics and quantum mechanics can be applied 
to the theory of chaotic spatially extended dynamical systems. An essential problem of the 
"physical" point of view is that one implies the nice invariant measure (or nice measures if the 
system is subjected to a phase transition) is defined. In general, coupled maps are not the case 
except the coupling parameter e is not vanishingly small, 0. 

Our example of coupled maps chain (|2|) has an invariant measure defined as a Gibbs state 
even if e is not small, fl7fl . The non-trivial collective behavior ("hidden coherence") is not 



actually captured by the present paper. At the same time, while computing the generalized 
entropy K q and the partition function Z, we shall perform the summations over " an 

configurations" {£}. The summation technique that we use relies upon the symmetry properties 
of the coupled chaotic maps chain and is known as a collective coordinates method. Proceeding 
with computations, we replace with an integration over an infinite set of continuous real 
valued positive variables. This stroke also implies that there is a number of nice measures on 
the space-time configurations. 

1.2 The scaling hypothesis for coupled map lattices 

The scaling hypothesis in thermodynamics has been an important milestone in formulation 
of the modern phase transitions theory. It is based on two fundamental principles of non- 
equilibrium thermodynamics: 

i) while interesting in critical phenomena, one can substitute a true micro-model for a fluc- 
tuation one; 

ii) in the regime of strongly developed fluctuations, the initial values are of no consequences 
for the system behavior. 

For the Ising model, the scaling hypothesis put forward by Kadanoff leans upon the 



"blocks" procedure, in which the interactions of spins have been replaced for the interactions 
of " spin blocks" . It consists of two propositions: 

j) the Hamiltonian of "spin blocks" interaction has just the same form as the Hamiltonian 
of spin interaction with parameters changed; 

jj) the dependence of parameters on the block size L has a scaling form as L — > oo. 

For the uncoupled maps having a "short range time memory", the scaling hypothesis has 
been propounded first in |ZD| and then used in [pT|| . 

We will follow a track by considering fluctuating effective quantities analogous to the Lia- 
punov exponents in finite time segments, n>l. The "effective Liapunov exponent" A™ define 



the rate of separation of points in tangent motion in finite time. They fluctuates around the 
true Liapunov exponent A having the fluctuation range A, [ 2"I |: 



n— >oo 



A = lim dA p„(A)A n , (6) 



in which p n (A) is the probability distribution to observe fluctuation A in n time steps. In the 
limit n — > oo, A n form a micro-canonical ensemble lim„^ 00 p n (A) = 6(A). 

An example of p n (A), one can consider the Fig.l where we have plotted the data of 5000 
numerical simulations for the uncoupled solitary map (H) at r = 1.6, for the effective Liapunov 
exponent A™, n = 700. The graph shows the probability distribution to observe deviations from 
the true Liapunov exponent A ~ 0.47 .... When n growths up, the bell-shaped profile becomes 
sharper. 

Being a product of n Jacobi matrices, the effective rates A n , for the maps with a short- 
range time memory, behave like averages of n random variables correlated only over short time. 
Therefore, their values, for n large, do not depend on the particular choice of initial conditions. 

The scaling hypothesis formulated in Jl9| for such maps claims on the self-similarity property 
for the microscopical entropy, S = logp(A) in a long however finite time string, n » 1, in the 
chaotic regime: 

p n (A)cxexp[-5], S = ng(A) (7) 

where g(A) is a "scaling function" (a "free energy") pertinent to the map, [ |H|j . 

In the Fig. 2, we have pictured out the ratio of scaling functions g n - L /g n2 v i & A computed 
numerically for the map (|j) in the sequences of n x = 600 and ri2 = 1200 time steps. The 
picture data show that there is an interval of deviations (~ [—0.1, 0.1]) within which the scaling 
function g(A) is almost universal at least in the first 1200 time steps. 

We note that in accordance with (0), given a smooth function F(A), the relevant observable 
is (F) = J dA p(A)F(A). Since n is a large number, one infers that the main contribution to 
this integral comes from the saddle-point configuration A* such that g{A*) is the maximal g. 

Indeed, this hypothesis breaks down for maps with strong long-time correlations, when a 
minor change of initial conditions can stir the behavior even in a very long time segments. 

In the present paper, we extend the scaling hypothesis to the case of coupled maps lat- 
tices. We search for the statistics of fluctuating quantities A£ analogous to the true Liapunov 
exponents A& in a long albeit finite time. We consider the extended lattices (N — > oo) of cou- 
pled maps with a short range time memory. By the way, the boundary conditions are of no 
consequences for the system behavior in "thermodynamic limit". 

Therefore, following [18] and resting on the numerical simulations data cited in [JI3], we 
propose that in the chaotic regime the "free energy" function gu{A) pertinent to the CML 
is universal up to the statistically independent block U, so that the probability to observe a 
deviation of range A = |A£ — A&| in the thermodynamic limit equals to 



p(A) oc exp[-S'], S = nN u gu{A). 



(8) 



Furthermore, we propose that in the chaotic regime the scaling function pertinent to the coupled 
maps lattice, gu, up to a change of parameters equals to that one for the solitary uncoupled 
map, g. 

In thermodynamic limit, according to (^|), the most important contribution to observables 
comes from saddle-points configurations A* (provided they exist) maximizing the scaling func- 
tion g u . 

1.3 Connections of a loosely coupled maps chain to a quantum mechanical system 
with a double-well potential 

Let us remark on a semblance between the statistics of effective Liapunov exponents of the 
simplest loosely coupled (e < 1) maps chain (|2|) in thermodynamic limit and of the quantum 
mechanical system with a double- well potential in the infinite (3 limit (i.e., in small temperature 
limit, T -> 0) |H. 

For |3 > 1, the main contribution to the partition function of the quantum mechanical 
problem comes from the instanton (or saddle-points) configurations which renders the action 
stationary. In the limit (3 — * oo, a linear combination of instanton solutions plays a role. The 
relevant contributions do not correspond to solutions of the classical equations of motion since 
these equations are non-linear. However, in the limit when all instantons are largely separated, 
such a combination also renders the action almost stationary since the relevant corrections 
coming from each instanton are exponentially small in the separation. 

To properly study the energy spectrum of the quantum mechanical problem at T — > 0, first, 
one has to introduce a set of collective variables (or coordinates) which describe fluctuations 
in distances between the instantons, and, second, to fix out their values in such a way that 
minimizes the variation of the action. 

In the present paper, we are going to apply this method to compute the asymptotic of a 
partition function Z over the space of possible configurations {c*.} generated by the chain of 
coupled maps (||). 

In Fig. 3.1, we have plotted the results of 10000 numerical simulations for the chain of 
N = 8 coupled maps (the value of coupling strength parameter has been taken as e = 0.02) 
supplied by the zero valued boundary conditions. This graph shows the probability distribution 
profile to observe the effective Liapunov exponents X n in the interval [0.44 . . . 0.48] inn = 1000 
iterations. One can see that almost all values A n drop around the true Liapunov exponents 
laying within this interval, A — 0.4700 . . . , Ai ~ 0.4660 . . . , and A 2 ~ 0.4548 . . .. The quota 
of A n arising apart from the close vicinities of Liapunov exponents A, decreases very fast as n 
grows up. 

In Fig. 3.2, we have drawn the distribution profile at n — 2500 time steps. For n ^> 1, 
the dispersion 5 becomes much smaller than the distances AAj between peaks, 5 <C AAj. This 
means that fluctuations which tend to change distances between Liapunov exponents induce 
just a vanishingly small variation of g at long time segments. 



We note that the thermodynamic limit when fluctuations of A£ are faded out, plays the role 
of the infinite separation limit for instantons in the quantum well problem mentioned above. 

The plan of the paper follows. In the Sec. 2, leaning on the thermodynamic formalism as 
applied to CML, we explain what kind of information on the Liapunov exponents spectrum one 
could obtain from the relevant partition function Z. Then, in Sec. 3 and 4, we compute Z using 
the scaling hypothesis formulated in the previous subsection (1.2) and following methods de- 
veloped in quantum mechanics with respect to the double- well potential problem |Z| . Namely, 
in the Sec. 3, we compute the scaling function g(A) (the free energy of the microscopical inter- 
action) for the map ([3]) coupled merely to the nearest neighbors. In the Sec. 4, we introduce 
the block scaling function <?w(A) (the block interaction) and then compute the asymptotic of Z 
in thermodynamic limit utilizing the collective variables method. We conclude in the Sec. 5. 



2 Thermodynamic formalism and the Liapunov exponents 



Thermodynamic formalism relies on the existence of a symbolic representation for the dynamics. 
The idea to study this representation via Gibbs states for the d + 1-dimensional system goes 
back to H and || . 

Let us suppose that one can split the phase space of CML LcZ d into disjoint boxes 6, of 
size a. We shall call as an orbit generated by the coupled map in site the sequence of boxes, 
{b\, . . . bj, . . . b n }, in which the coupled map takes values at consequent time steps j = 1 . . . n. 
We shall call as a configuration the set of orbits, 

Cl = {{h,...b n } Xi ,Xi G L} , 

in the finite piece L C Z d . The thermodynamic formalism comes about by asking for a con- 
ditional probability P f ££| ^° observe the particular configuration given the config- 
uration in the complement £2 d \z/ This probability looks like an expectation in the canonical 
ensemble, 

(9) 



exp 



*Z d \L 



P{C L \C 

in which the normalization factor (the partition function, in terms borrowed from statistical 
physics) is 

m 

(the sum is taken over all configurations depending upon possible configurations Q d \ L in the 
complement). Here, f3 > is an inverse temperature parameter representing a "heat bath" of 
coupling maps. 

The interaction potential $ plays the role of a Hamiltonian defined for the finite sublattice 
Lby 



1L 



For infinitely long time strings, it depends upon the configuration £ Z d in the whole lattice. The 
Jacobi matrix (DF) can be found in the form (DF) = B(I—W) where B is the diagonal part of 
the matrix which provides a point-wise contribution from maps with no coupling, and W comes 
from coupling. The coupling W possesses the translation invariance property. Therefore, the 
potential ( pTTj) can be extended naturally to all subsets of the lattice Z d by translation invariance. 

For future convenience, we use the identity det = exp Tr log to rewrite the potential ([TT[) in 
the following form 



log(£F)(££|£^)], (12) 



which has been discussed extensively in ||-[25|, and [ 26| . Due to expansion log(DF) = logB 



log(J — W), the Hamiltonian flT2"| ) can be divided in two folds: $£ = $o + U, where $o is a 
"free" Hamiltonian relevant to the uncoupled maps and the coupling "interaction" U. 

Using the definition of effective Liapunov exponents, A£, one can express the partition func- 
tion (pi]) as follows, 

Z = lim V exp (nPTr[X£\) = exp ((3Tr[X k }) (13) 

n— >oo ' J 

where the sum * s taken over the ensemble of fluctuating quantities A'£. For the uncoupled 

maps, the Eq. ( |T3| ) yields A = lim^oo log Z/nfi. However, for coupled map lattices, one can 
just obtain merely the sum of positive Liapunov exponents Tr[AJ. 

It is essential to make use of the scaling hypothesis (Sec. 1.2) while searching for the 
Liapunov exponents. In particular, it implies (see Sec. 4) that 

Tr[A fc ] ~ N>1 N. (14) 

Therefore, at least for large N, the values of the Liapunov exponents Xn can be found by the 
simple subtraction, A^v = Tr[AAr_i] — Tr[Ajy]. Furthermore, the maximal value in the spectrum, 
Ao equals to the Liapunov exponent of the uncoupled map and, at least for the chain (|), is 
also known. 

However, since the relation (|I4D is formally valid only for N ^> 1, the spectral values 
Afc, k ~ 0(1), are still uncertain. The numerical simulations of systems demonstrating the 
extensive chaotic behavior (i.e., scaling) give an evidence (see for example, [^]) in favour of that 
the threshold lattice size is typically N ~ 8. 

3 Scaling function g(A) for the map coupled merely to the nearest 
neighbors 

The scaling function g(A) plays a role of free energy in the ensemble of fluctuating effective 



Liapunov exponents. In the limit n — ► oo, it can be related to the generalized entropies [19 
K q defined by 



K q (A) = -J- lim ±log£> n (A)*, (15) 
q — 1 rwoo n L — ' 



via Legendre transforms, [19 

1 



KM 



T [Aq-g(A)], g(A) = X (q) 



9 ( \ 
Q—X{Q), 



X (q) = {q- l)K q . 



(16) 



q — I ~ ' ' " " dq 

It has been extensively discussed in a literature that lim ? ^ K q and lim ? _,.i K q are the topological 
and metric entropies, respectively. An example of computation of g{A) has been given in W§ 
for the simple uncoupled two-dimensional "baker's" map. 

Let us choose the site as an origin and introduce a point-wise Hamiltonian 



x E0 + m(o,v o ), 



(17) 



in which the interaction potential U^(0, Vq), and Vq is a nearest neighborhood of 0. In the 
present section, we suppose that there is no interactions between Vq and remainder of the 
lattice. 

The free point-wise Hamiltonian equals to the logarithm of diagonal elements of the Jacobi 
matrix, A e0 = log (r(l — 2e)) . The potential U^(0, V ) describes the instantaneous fluctuations 
of the effective Liapunov exponent. It comprises of a sum of contributions coming from the 
neighboring sites y e Vo, 

c^(o,y ) = X)<o(y). 

yev 

The relevant point-wise generalized entropy takes the form: 



K q (w e 



lim — log exp 

— 1 n-HX n ^— ' 



\w £ - f3q <o(y) 
yeVb 



in which w £ = Ao — A e o, and X]{A"} denotes the sum over all possible orbits {A n }. Note, that 
the first term in the exponent, fi = exp [(3qw £ \ persists for any orbit. It plays a role of the 
"fugacity" of the fluctuation "gas". 

To perform the summation in properly, we note that the site and its nearest neighbors 
Vo form an isolated system. Therefore, any fluctuation risen in can be represented as a 
vector in the V-dimensional space of particular contributions from the neighboring sites, U = 
{u e n o(y)}. Since the entropy of closed system does not change in time, the sum ^{A"} can 
be interpreted as an integral over the set of linear transformations U — ► U which render the 
entropy fll8D unchanged. 

The transformations we are interested in is apparently defined by the set of positive param- 
eters ify > 0, such that 

Vn. (19) 

yeV 

Choosing at each time step V — 1 values {(fy} at random, one generates a particular orbit {A n }. 
Therefore, one can replace the sum over possible orbits in flT8p with integrations over the real 
positive variables ip™, 



K q {w £ ) 



J ^v >0 n&+ ydV Q 



X 



(20) 



x e n g [/3roe-l/nE ye v ^<o(y)] , § | p _ 



yev 

The factor 1/n before the integral over ip™ in ( p0|) arises because the result of integration has 
to be invariant under a cyclic permutation of ^"—variables in time. 
After the integration over ip™ variables ||27|| ) in (0), one obtains 

-rj+ioo e P(s+nvj £ ) 

lim — log > / 

q 



K q (w £ ) 



lim — log / 

— 1 n^oo n Z — ' 2Tcni I _ 



-cfe, 



(21) 



in which (s,q) = Yl n( zz + Ylyev ( s ~ q u no(y)) is a the polynomial of the order V. For 
piecewise linear maps, one can proceed further since the slope of map is a piecewise constant, 
and the values u n0 (y) in the polynomial Py\s, q) does not depend on n. 
The sum J2nez+ therefore can be easily performed, 

-r/+ioo 

cfe/3e _/3s log w(s, zu £ 



K q {w £ ) 



1 



q-1 



log 



2m(3 



77—400 



22 



where we have defined 



Wis, W £ ) = 1 



The remaining integration in (|22| ) is performed by parts, 



K q {w £ ) 



q-1 



log 



2m(3 



v -ioo W{S,W £ ) 



The contour can be deformed in the half-plane dt(s) < to enclose the poles of the integrand 
given by the solutions of the equation 



II (8-<p% (y))=e-**. 
yeVb 

The solutions of the equation ( p3| ) can be easily found, 

(3 



(23) 



s = qu £ n0 + exp[-— w E ) 



Then, the residue theorem yields: 

K q {w e ) 



q-1 



q— + qu n0 - e 



Using the Legendre transforms (|TBp, one obtains the relevant scaling function, 



pe- 



rn 



(25) 



(26) 



The scaling function g(A) (BB) has a typical form pertinent to one-instanton contributions as 
they appear in different problems in quantum mechanics. We also note the relation between 
the fugacity of fluctuation gas /i and g(A), fi(A) = g v (A). 



4 The partition function in the thermodynamic limit 



In accordance to the scaling hypothesis, while proceeding to the interactions between blocks, the 
Hamiltonian of the interaction has to be the same as the microscopic one, except the values of 
parameters. The general scaling in the system (i.e., the fact that the dependence of parameters 
on the block size L has a scaling form as L — > oo) then follows as a natural consequence from 
the latter preposition p8| . 

Proceeding from (|26|), for the statistically independent block U (which constructive definition 
we give below), one obtains the relevant interaction, 



U 



(27) 



9u(A) = exp 

It differs from (^) just on the volume U. 

Let us consider a finite piece L of the lattice 1 d . We assume that L is large enough to 
comprise of many statistically independent blocks U. For each site x G L, we assign an inverse 
temperature parameter 9 X > 0. These parameters is to be subjected to two obvious constrains 
expressing the idea of statistical independency of the block from the rest of lattice: 

1. The statistically independent block has to be in thermodynamical equilibrium with the 
remainder of lattice. Therefore, 

N U 

= S>. (28) 

5=1 

where the summation is going over the sites belonging to the same block U. 

2. All parts of the sublattice L should be in thermodynamic equilibrium with respect to each 
other, 

u = w u p*-' < 29 > 

in which U denotes the number of lattice sites belonging to one statistically indepen- 
dent block, iV is a total number of sites in L, and Nu is a total number of statistically 
independent blocks in L. 

Pursuant to (^), the block- wise partition function has a scaling form, 

Z = ^exp[-/3A/ W ( 7w (A)] (30) 

{ u} 

where the summation is performed over all statistically independent (in the sense of constrains 
1.) and 2.) ) collections of sites {U}. To compute this sum in the thermodynamic limit 
(Nu — > oo), we employ the collective variables method again. 

In terms of the ^-variables, the block-wise scaling function reads as 

N 



gu(A) = exp 



-—ye. 



and therefore, 
Z 



Nu 



E 

k=l 



e Xa >o 



H d6 * 



exp 



-(3N u e N u 



Nu 



s=l 



(31) 



where the factor 1/k in front of the integral is risen due to cyclic permutation of variables. 
We use the integral representation for the 5— function ( ||27|| ) and introduce the new auxiliary 
variable a such that 



Nu 



with the definition: 



E 

k=l L 



J (a) 



1 

2~i 



Nu 



-ri+ioo J '« i 

da e^^-r [J(a)] A 



—rj—ioo 



k=l 



(32) 



exp[a6 + (3N u e eA/Nu ] d9. 



By giving to a a small negative real part, we have insured the convergence of the integral J {a). 
By means of the change of variables: r = /3Nue~ 9m£ / Nu , one obtains, 



J(a) 



pNu dr f(3N u 

T 



In thermodynamic limit Nu —> oo, this integral converges to the T— function, 

J(a) ~ («*r (^) . 



The sum over k in (|32| ) can be performed easily: 

— / da log 



1-{(3N U )^T 



N u a 
A 



(33) 



The integral in (|33|) can be computed by parts. Due to asymptotic behavior of the T— function, 
the integral converges, and the contour of integration can be deformed to enclose the poles 
of the integrand in the half-plane Re(a) < 0. These poles correspond to the solutions of the 
equation: 



((3Ni 



r 



Nua 
A 



(34) 



In the limit Nu — > oo, solutions of the equation fl3~4]) are close to the poles of T— function, i.e. 
to negative integers. 

Since the coupling strength parameter e is small, one can taylor the T— function in powers 
of e to find out the approximate solutions of (p4|). In the first order in e, this procedure gives: 

oliNu 



-Is 



1...N. 



Finally, one arrives at the spectrum of possible values 

2e£ 



N 



u 



The asymptotic result for the spectrum of Liapunov exponents follows from (|35|), 



2s 



1...N. 



(35) 



(36) 



The spectrum (|36|) is consistent with the relation (|5[) observed numerically for wide variety of 
CML. One can see that for ^> 1 and e small it resembles closely the spectrum (01). 



5 Discussion and Conclusions 



In the present paper, we have proposed and analyzed the statistical theory of fluctuations of 
the effective Liapunov exponents which correspond to the true Liapunov exponents in finite 
time segments. 

We have formulated the scaling hypothesis for the CML and justified the general scaling 
behavior (|5|) of the CML in the extensive chaotic regime. 

Due to presence of large parameters (n and N), fluctuations of the effective Liapunov expo- 
nents can be well described by the saddle-point approximation. 

In the framework of thermodynamical formalism in CML, we have computed the scaling 
function g(A) for the 1— dimensional chain of coupled piece-wise linear maps and demonstrated 
that it has a form typical for the 1— instanton contribution in the quantum mechanical problem 
with a double-well potential. 

The basic technical point of scaling hypothesis is that the block-wise Hamiltonian appears 
to have the same form as the micro-interaction. As a direct consequence of this assumption, 
one arrives at the following scaling properties of the Liapunov exponents spectrum for large 
lattices, iV > 1: 

N\ k>0 ~ N/U, and Tr(A fc ) ~ N. 

They allow to compute the asymptote of the Liapunov exponents spectrum for N ^> 1. 

We should note that the developed technique can be also applied to the multidimensional 
coupled piecewise linear maps. However, it probably fails in case of non-linear map lattices. 
Nevertheless, since most of non-linear maps can be approached by piecewise linear ones, we 
believe that our technique would also give good approximations to more complicated loosely 
coupled non-linear map lattices. 
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CAPTION FOR FIGURES 



Figure 1. 

The probability to observe deviations of A n , n = 700, from the true Liapunov exponent 
A w 0.47. . . (5000 simulations). 
Figure 2. 

The ratio of scaling functions g ni /g n2 Yia - the deviation range A computed in the sequences 
of rii = 600 and n 2 = 1200 time steps. 
Figure 3.1 

The probability distribution profile to observe the effective Liapunov exponents A™ in the 
interval [0.36 . . . 0.46] in n = 1000 iterations (10 4 simulations, the chain of N = 8 coupled maps, 
the coupling strength parameter, e = 0.02). 

Figure 3.2 

The probability distribution profile to observe the effective Liapunov exponents A™ in the 
interval [0.36 . . . 0.46] in n = 2500 iterations (10 4 simulations, the chain of N = 8 coupled maps, 
the coupling strength parameter, e = 0.02). 



FIGURE 1 



0.3H 

c 
c 



S 0.2 
30.15 



00 





v 
o 



0.4 0.42 0.44 0.46 0.48 0.5 0.52 0.54 
Effective Liapunov exponents 



0.2 



0.15 





1 





1 







































o 

* 


: V 













# 




0* 





H 



0.455 0.46 0.465 0.47 
Effective Liapunov exponents 



1.8 



s 1.6 



(9 " 

(9 - 

'5 - 

I . 

\ 1.2 



FIGURE 2 









1-| ^00000000000^0^0 



41 405 0.05 0.1 0.15 0.2 



0.25 




FIGURE 3.2 








% 




0.2 






<o 
<b 
<fc 




0.15 






00 
00 
00 




0.1 


«> 




00 






oo 













00 


00 




0.05 


o« 


00 




00 

















! ( > 


00 


8 







o H 




n o 






0.455 _ 


0.46 0.465 


0.47 




Eff 


ective Liapunovexpone 


its 



